ID Exciton Spectroscopy of Semiconductor Nanorods 
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We have theoretically shown that optical properties of semiconductor nanorods are controlled 
by ID excitons. The theory, which takes into account anisotropy of spacial and dielectric confine- 
ment, describes size dependence of interband optical transitions, exciton binding energies. We have 
demonstrated that the fine structure of the ground exciton state explains the linear polarization of 
photoluminescence. Our results are in good agreement with the measurements in CdSe nanorods. 

PACS numbers: 73.22.-f, 78.67.-n, 71.70.Gm, 77.22.Ej 
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There is growing interest in nano-size crystalline semi- 
conductor structures of various shapes such as nanocrys- 
tals (NCs) nanor ods (NRs) @ and nanowires (NWs) 
Q created by the "from the bottom up" technological 
approach. Size-tunable control of their optical and trans- 
port properties combined with the ability to move these 
nano-size objects around with precise control opens the 
exciting possibilities for the creation of new functional 
materials which can be used in unlimited applications. 
Among these nanostructures, the NCs are the most heav- 
ily studied and one can find a broad description of their 
properties and their potential applications in the reviews 
of Brus0 and Alivisatos0|- 

The optical properties of NRs, however, differ signifi- 
cantly from those of NCs. Compared to NCs, the NRs 
show higher photoluminescence (PL) quantum efficiency 
0, strongly linear polarized PL0,0|, an increase of the 
global Stokes shift 6] , and significantly faster carrier re- 
laxation |8|]. The Auger processes in NRs are strongly 
suppressed relative to those in NCs 0, 0] , which subse- 
quently decreases the optical pumping threshold for stim- 
ulated emission [lllilt. The size and shape dependence 
of the optical and tunneling gaps measured in CdSe NRs 
(l2^ shows an unexpectedly large difference that cannot 
be explained by the electron-hole Coulomb correction to 
the optical gap used for NCs. To describe their measure- 
ments Katz et al. applied the four-band theory devel- 
oped by Sercel and Vahala The most important re- 
sult of the latter paper was the prediction of an inverse or- 
der of light and heavy hole subbands in NWs, which was 
later confirmed by semiempirical pseudopotential calcu- 
lations 0, llil flfl and numerical calculations within the 
6 band model [lq . An anisotropy of the light hole-to- 
electron optical transition matrix element predicts some 
deg ree of the PL linear polarization parallel to the NW 

The dielectric confinement connected with the differ- 
ence between dielectric constants of semiconductor crys- 
tallites, k s , and surrounding medium, K m , (for example, 
see [l9() usually does not affect the optical spectra of 
NCs because the charge distributions of an excited elec- 
tron and hole practically compensate each other at each 
point of the NC. This leads to a complete screening of an 



electric field of the total charge from penetration into the 
surrounding medium. This is not the case, however, in 
NRs where the electron and hole are at a distance larger 
than a NR radius and interact predominantly through the 
surrounding medium, which usually has a small dielectric 
constant, K m <C k s . This should lead to a formation of 
ID excitons (IDEs) with large binding energy similar to 
that in NWs [iO] . In addition, the dielectric confinement 
of anisotropic crystallites leads to strong PL polarization 
and PL polarization memory that was observed in NWs 
[2l| and in porous Si j^- This dielectric model, how- 
ever, does not explain the strong linear polarization of 
NRs @. 

In this letter, we have developed a theory that de- 
scribes energy spectra and polarization properties of 
IDEs in NRs. Our theory takes into account both spa- 
tial and dielectric confinement and allows one to describe 
optical properties of elongated NRs with quasi-cubic and 
zinc-blende lattice structure. The calculation for CdSe 
NRs surrounded by a dielectric with n m = 2.0 has shown 
that the IDEs control optical properties of NRs even at 
room temperatures and has described the difference be- 
tween the optical and tunneling gaps 01 • We have shown 
that the fine structure of IDE explains the observed 87% 
polarization degree of PL |(J and leads to the fast radia- 
tive decay in NRs at room temperatures. 

In what follows, we consider the NR as a crystalline 
that has the shape of an ellipsoid of revolution with ma- 
jor semi-axis b (the NR axis) much larger than the minor 
semi-axis a (the NR radius). The motion of electrons 
and holes is analyzed within the 6-band model that has 
successfully described optical spectra of CdSe NCs |23| . 
The significant elongation of NRs (b 3> a) makes it possi- 
ble, by using an adiabatic approximation, to separate the 
parallel motion from the motion perpendicular to the NR 
axis. Within this approximation we initially neglect the 
slow motion of carriers parallel to the NR axis. Assum- 
ing that the momentum of the parallel motion k z ~ 0, 
we find the spectrum of carriers in 2D confinement per- 
pendicular to the NR axis. Next, we consider the parallel 
motion by averaging the Hamiltonian over the fast mo- 
tion of carriers strongly confined in 2D. 

Each electron state in NRs is characterized by the 
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angular momentum projection on the NR axis, m = 
0, ±1,.... The wave functions of the electron motion 
perpendicular to the NR axis can be written as ~ 
e %m ^ J| m | (kp), where (j> is the azimuthal angle, p is the 
distance from the NR axis, and Ji(x) is the Bcssel func- 
tion of the Z-th order. The momentum hk is connected 
to the kinetic energy of the electron motion perpendic- 
ular to the NR, E e = h 2 k 2 /2m(E e ), where m(E) = 
mo/{l + 2/ + (E p /3)[2/(E + E g ) + l/(E + E g + A)]} is 
the electron effective mass that depends on the electron 
energy, E, mo is the free electron mass, E g is the band 
gap, A is the spin-orbit splitting of the valence band, E p 
is the Kane energy parameter, and the parameter / in- 
cludes the contribution of remote bands. Assuming that 
the wave function vanishes at the NR surface, we find 
the momentum Tik which satisfies the boundary condi- 
tion: ^fc„| m | = ha n \ m \/a, where a„| m | is the n-th root 
of the Bessel function of the |m|-th order. This gives 
the following energy of electron subbands E n \ m \ with the 

corresponding wave functions "i"^, (p, </>): 
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where p — p/a and the three smallest a„| m | are a±o ~ 
2.41, an « 3.83, a 12 « 5.14. 

Due to a strong spin-orbit coupling in the valence 
band, each hole state in the NR is characterized by 
the total angular momentum projection on the NR axis 
j z = m + J z which is the sum of angular momentum 
projection and the projection of the hole spin J = 3/2 
(j z = ±1/2, ±3/2, ...). The energy spectrum of holes is 
found within the 6-band model, which takes into account 
the nonparabolicity of light hole spectrum j^]. In this 
model, the hole wave function of the motion, perpen- 
dicular to the NR axis, with total angular momentum 
projection, j 2 , can be written 

M =-3/2 V 2 * 



1/2 

E 

m=-i/2 



E C J' tj \J,-^p)^=^u^ (2) 



l % ^ and ^ are the valence band Bloch func- 



where u\ / 9 

tions |24|. For every value of the hole energy E, there 
are dimensionless momenta Xt and the sets of coefficients 
C^ : * corresponding to the three branches in the disper- 
sion law of the valence band (t = I, h, s): 
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where e = E/E , 5 = A/E , E = fr 2 /2m a 2 , and the 
energy dependent Luttinger parameters 7 f and 7 L are 
connected with the contributions from remote bands, 7 i 
and 7 , by the relationships: 7 f =_7i + E p /[3(E g + E)\ 
and 7 L = 7 + E p /[6(E g + E)] 24]. At k z = 0, the 
6-band Hamiltonian can be decomposed into two 3- 
band Hamiltonians and the wave functions 5 , '' ± (p) in 
Eq.(2) are described by the set of only three coefficients: 
C u± ee (C^y^C^y'lC^'l). For the heavy hole and 
for the mixture of the light- and spin-orbit split holes, 
these coefficients are C K± = (\/3,l,0) and C sW,± = 
(l,-V3,Xa(l)/iV$) respectively, with Xs(i) = 1 - 
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FIG. 1: Size dependence of the electron and hole energy spec- 
trum in CdSe NRs [E g = 1.839 eV, E p = 19 eV, A = 0.42 eV, 
/ = —1.035, 7 = 0.55, and 7 i = 2.1). Arrows show optically 
allowed interband transitions. Inset: Size dependence of the 
hole effective mass at the l£i/2 subband. 

The hole wave function of the motion perpendicular 
to the NR axis is a linear superposition of the three 
eigenfunctions ^f 1 , ^fjS an d ^j'^' corresponding to 
the same energy (see Eq.(3)). The requirement that the 
wave function vanishes at the NR surface defines the coef- 
ficients of this linear superposition and gives the disper- 
sion equation for the hole quantum size levels (QSLs). 
The wave function can be written: 
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l)/[V3J\ jzTl /2\(x h )}, 
and the coeffi- 



3zh - tf |j»Tl/2| 
U s -XlJ\j^l/2\(xi)/[XsJ\j zT l/2\(Xs)}, 

cient is determined by the normalization condition: 
/ d 2 p\^f z (p, <t>)\ 2 = 1. The energy of the hole QSLs are 
given by the solutions of: 

3#- (x a ,x h ,xi) - 3g i (xi,Xh,x B )+g i (xh,x 3 ,Xh) 



gf(x h ,xi,x h ) = 



where gj^ («, v, w) = v 2 [u 2 — (7^ + 27^ 
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The size dependence of the electron and hole QSLs 
calculated for CdSe NRs is shown in Fig.l, where we 
use level notations adopted from molecular physics. For 
the electron levels it is nA e , where A = S, II, A, ... for the 
states with angular momentum projection on the NR axis 
\m\ = 0, 1, 2, .... For the hole it is nAu i, where A corre- 
sponds to the smaller value of \m\ and |m + 2|, the abso- 
lute values of the two angular momentum projections rep- 
resented in the hole wave function, and n is the number of 
the level for the given symmetry. The 6 band model cal- 
culation shows that the IS1/2 hole level with \j z \ = 1/2 
is above the IS3/2 level with \j z \ = 3/2 in consistence 
with previous predictions of the reverse order of the light 
and heavy hole QSLs in NRs. The arrows in Fig.l show 
the optically allowed transitions between the electron and 
hole QSLs, which have been obtained by the calculation 
of the overlap integral K nA>n , A >^ = J d 2 p^{ (p)^(p) 
(see Ref.[23l|L The energy shift 8e nAj of each ID hole 
subband at finite k z can be described in terms of the ID 
effective mass: 5e nAj (k z ) — % k 2 /2m nAj . The depen- 
dence of the first hole subband effective mass, mis 1/2 , on 
the NR radius is shown in the inset of Fig. 1 . 

The calculated optical spectra, however, does not in- 
clude the electron-hole (e-h) Coulomb interaction en- 
hanced by the penetration of an electric field into a sur- 
rounding medium. These e-h interaction can be written 

U(r e ,r h ) = -e 2 / (k s \r e - r h \) - eV s (r e , r h ) 

+ eV s (r e , r e ) /2 + eV s (r h , r h ) /2 , (6) 

where the first term is the direct e-h interaction in a semi- 
conductor, the second term is the e-h interaction medi- 
ated by the surface separating the NR from the surround- 
ing dielectric, and the last two terms are self-interactions 
of each particle with the surface. For an ellipsoid, the 
potential V s can be written in the following form |25j 
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X K l W) Pn (V) Pn (£') Pn (0 COS TU - <j>) , (7) 

where P™, Q™, P™, and Q™ are Legendre functions |2^| 
of spheroidal coordinates (77 and £) and their derivatives, 
Sij is the Kronecker symbol, and 
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The spheroidal coordinates are connected with cylin- 
drical coordinates (z and p) by the relationships: 

irt^/W^ = z and [(1 - £ 2 ) ( V 2 - 1) (b 2 - a 2 )] 1/2 = p. 
At the spheroidal surface of the NR 77 becomes r/ s — 
b/Vb 2 - a 2 . 

In the adiabatic approach, one averages the potential 
U (r e , rh) in Eq.6 over the electron and hole wave func- 
tions of the fast motion perpendicular to the NR axis 



(Eq.l and Eq.4). The procedure leads to the ID potential 
of e-h Coulomb interaction and the shift of the electron 
and hole QSLs due to their interaction with the "mirror 
forces" described by the last two terms in Eq.6. The ID 
potential t/is c -is 1/2 i z ) for the lowest lS e electron and 
l£i/2 hole subbands can be approximated as 

UxE.-XE 1/a (z) ~ C4ff M = ~e 2 /[n m (\z\ + p cS )} (9) 

where p s — 0.7a for n s — 6.1. This approximation is 
valid only if the distance \z\ = \z e — zh\ between two 
charges is much smaller than the NR length 26. 

In the ID adiabatic potential, IDEs are formed under 
each pair of e-h subbands. These IDEs strongly modify 
the absorption and PL spectra of NRs. The approxima- 
tion of the ID potential in Eq. 9 makes it possible to ana- 
lytically find energy spectra and wave functions of IDEs. 
The wave function of the IDE is 4> ex (z) = AW a ^/ 2 (z), 
where W a ,i/2 is the Whittaker function [2(|, A is the 
normalization constant, and z = ±2 (\z\ + p c g) / (aica) 
with a\£> — e 2 K m /h 2 p, where p — [l/m(E) + l/m n \ ^]~ 1 
is the reduced electron and hole effective mass. The plus 
sign is for positive z and the minus sign is for negative z. 
Each state of the IDE is characterized by its parity. For 
"even" states, including the ground state, the derivative 
of the wave function must turn to zero at z = 0. The wave 
function of "odd" states becomes zero itself at z = 0. The 
boundary conditions determine the n-th value of a n and 
the energy spectrum of IDEs as E n — — h 2 /(2paf D a 2 l ). 

The size dependence of the CdSe NR transport gap, 
which is determined by the energy between the lowest 
lE e electron and IS1/2 hole subbands, is shown in Fig. 2 
for the bare QSLs (dotted line) and for the QSLs shifted 
by the "mirror forces" (dashed line). One can see that 
the "mirror forces" corrections to the transport gap of 
NRs are significant. The optical gap in NRs is deter- 
mined by the IDEs. The size dependence of the first two 
"even" IDE states (only "even" states are optically ac- 
tive) is shown in Fig. 2 by solid lines. Even in thick NRs, 
the binding energy of the ground and excited states is 
higher than 150 meV and 50 meV, respectively. As a 
result, IDEs should completely control the absorption 
and PL of CdSe NRs even at room temperature. Given 
the good agreement with the experimental data for the 
transport and optical gaps in CdSe NRs [l^, formation 
of the IDEs in NRs explains the substantial difference 
(~ 240 — 300 meV) between the two gaps in NRs with 
a = 1.8 - 2.0 mi. 

The above results for IDEs in NRs can be applied 
directly to NWs. Unlike in NWs, however, the ra- 
dius of NRs decreases and causes a rise of the elec- 
tron and hole QSLs energy. Both electrons and holes 
tend to be at those parts of the NR where the radius 
is larger. In NRs of ideal ellipsoidal shape, this leads 
to the parabolic potential, acting on the IDE center of 
mass motion, whose minimum is at the NR center and 
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FIG. 2: Size dependence of the transport and optical energy 
gap in CdSe NRs. Dotted and dashed lines show bare and 
dielectric confinement dressed energy between the IE 1/2 hole 
and lE e electron subband. Two solid lines show the optical 
transition energy for the two first optically active IDEs. Ex- 
perimental data for the transport and optical gap measured in 
CdSe NRs from Ref. |T3| are shown by empty circles and filled 
squares. Insets a: Dependence of the IDE localization energy, 
HQ/2, on the NR semi-axis length b. b: Size dependence of 
the ground IDE fine structure. 



frequency is O = Ja[—d{Exo + E-^ 1/2 )/da]/b 2 M 1 where 

M = m(-E) + mi£ 1/2 is the IDE mass. The inset in Fig. 2 
shows how the IDE localization energy, Ml/ 2, depends 
on b for two values of the NR radius. 

The PL properties of NRs are determined by the fine 
structure of the ground exciton state. The electron-hole 
exchange interaction H cx = — (2ir /3)hujsT(a<ex)' i fi( r e — 
r h )(s-J) [53, where s is the electron spin 1/2 matrix, a cx 
is the bulk exciton Bohr radius, and Kust is the bulk ex- 
citon singlet-triplet splitting (in CdSe: Tilost — 0.13meV 
and a cx = 56 A), splits the four- fold degenerate ground 
state of the IDE in NRs into three states: two states with 
total angular momentum projection F z = s z + j z = 
and a degenerate state with F z = ±1. The calculation 
of the IDE fine structure is similar to Ref. (2^1 and it 
shows that the ground d IDE is the optically passive 
state with F z — 0. The other two IDE states, ±1 with 
F z = ±1 and b with F z = 0, are optically active. The 
size dependence of this fine structure is described as: 

e d j0 b = e ex [r)[f3) T \(0)], £±1 = ~£ C yJl{P) (10) 

where £ cx = (tt / '3) fiuj S t (al x /a 2 ) $ 2 X (0) and r)(f3), A(/3) 
are the dimensionless functions of the ratio of the light 
and heavy hole effective masses (3. In CdSe NR's r\ {(3) = 
0.34 and A(/3) = 1.38. The dependence of the exciton 
fine structure on a, is shown in the inset (b) of Fig. 2. 

For the exciton localized at the NR center, we 
have calculated the dipole matrix element of the op- 
tical transition to the b exciton state: d b = 
v /2E p /3m uj 2 K 1SA s 1/2 (irh/4:Mn) 1 / i 'f> cx (0), where u is 
the optical transition frequency. This dipole is 2^/2 times 



larger than the d±i dipoles of the optical transitions to 
the ±1 exciton states. In the ideal NR, the exciton PL 
intensity is proportional to the sum of the relative popu- 
lation of the ±1 and b exciton states multiplied by the 
corresponding matrix element squared. In addition, due 
to the NR anisotropy, the emission from the b state is 
much stronger than from the ±1 states because the elec- 
tric field of a photon is significantly reduced if the field 
is perpendicular to the NR axis while the photon field 
remains almost unchanged if it is parallel to the NR axis. 
This dielectric enhancement can be written as 

R c = {[K m + (k s - K m )rc (X) ]/[ K m + («„ - Km)n (Z ']} 2 , 

(11) 

where = a 2 [ln(2&/«) - l]/b 2 and = (1 - n^)/2 
are the depolarization coefficients for the two components 
of an electric field parallel and perpendicular to the NR 
respectively. In the elongated (6 ^> a) CdSe NRs R c w 
)/2K m ] 2 « 4.1. Taking into account the exciton 
level population, the anisotropy of the matrix elements, 
and the dielectric confinement, we have calculated the PL 
polarization degree of a single NR as P = (7ii — (/y + 
where 7m and I± are the intensity of light with the 
polarization vector lying in the same plane with the NR 
axis and perpendicular to it jy| . The angular dependence 
of P can be written: 

P = cos2^[(l-x)sin 2 6']/[sin 2 6' + .x(l + cos 2 6i)] , (12) 

where x = e^ e o i > _e ± 1 )/' cT /(4i? c ), T is the temperature, ip 
is the angle between the polarization analyzer and the NR 
axis, the observation angle 9 is measured from a direction 
parallel to the NR. At room temperature and 8 = ir/2, 
the polarization is modulated at 87%, which is in very 
good agreement with the experiment |fj. 

Our calculations show that the anisotropy of dielectric 
and spatial confinement shortens the radiative lifetime at 
room temperatures down to 500 ps in reasonably elon- 
gated NRs relative to ~ 20 ns measured in CdSe NCs. 

The IDE is localized at the NR center only if the 
NR thickness decreases to the edges. Any nonmonotonic 
variations of NR thickness, NR bending, and localized 
charges at the NR surface should lead to IDE localiza- 
tion at different parts of the NR. This results in inhomo- 
geneous broadening of IDEs, even in a single NR. 

To summarize, we have described the structure of the 
ID excitons in NRs, which control practically all aspects 
of the NR PL properties even at room temperatures. 
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